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Abstract. We present a new approach to analyze the validation of weakly nonlinear geometric 
optics for entropy solutions of nonlinear hyperbolic systems of conservation laws whose eigen- 
values are allowed to have constant multiplicity and corresponding characteristic fields to be 
linearly degenerate. The approach is based on our careful construction of more accurate auxil- 
iary approximation to weakly nonlinear geometric optics, the properties of wave front-tracking 
approximate solutions, the behavior of solutions to the approximate asymptotic equations, and 
the standard semigroup estimates. To illustrate this approach more clearly, we focus first on the 
Cauchy problem for the hyperbolic systems with compact support initial data of small bounded 
variation and establish that the L 1 — estimate between the entropy solution and the geometric 
optics expansion function is bounded by 0(e 2 ), independent of the time variable. This implies 
that the simpler geometric optics expansion functions can be employed to study the behavior 
of general entropy solutions to hyperbolic systems of conservation laws. Finally, we extend the 
results to the case with non-compact support initial data of bounded variation. 



1. Introduction 

We are concerned with weakly nonlinear geometric optics for entropy solutions of the following 
hyperbolic system of conservation laws: 

(d t U + d x F(U) = 0, UGW 1 , 

\U\ t= o = U° + eU 1 (x), { - } 

where F: W 1 — > M. n is a smooth function. The Jacobian matrix X7F = VuF{U) is diagonalizable 
at any point U £ K n and has n real eigenvalues such that any eigenvalue has constant multiplicity. 
Without loss of generality, we assume that 

Xi(U) < ■ ■ ■ < \ m (U) < \ m+ i(U) = ■■■ = X m+p (U) = X(U) < X m+p+1 (U) < ■■■< X n (U), 

(1.2) 

where 1 < m < n — 1, 1 < p < n — m, and the corresponding left and right eigenvectors 
{1, ([/)}? =1 and {r k {U)} n k=1 satisfy 

lj(U) VF(U) = XjljiU), VF(U) r k (U) = X k (U)r k (U), \,{U) ■ r k (U) = S jk . (1.3) 

Here dj k is the Kronecker delta satisfying 6jj = 1 and Sjk = when j ^ k. When p = 1, the 
system in (jl.ip is strictly hyperbolic. In addition, we assume that each characteristic field is 
either genuinely nonlinear or linearly degenerate, that is, either of the following holds: 

VuXj(U) ■ Yj{U) = 1 or VuXj(U) • tj(U) = for all U G R n , j = !,■■■ , n. (1.4) 
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A prototype of such hyperbolic systems is the full steady Euler equations for supersonic ideal 
gases in M 2 with the following form for x£K 2 : 

V x • (pv) = 0, 
V x • (pv <g> v) + Vp = 0, 
V x • ((ip|v| 2 + ^L)v) = 0, 

where p is the density, v = (i>i,i>2) the fluid velocity, and p the pressure. The notation a ® b 
denotes the tensor product of the vectors a and b. 

Set U = (pv\, pv iv, /wi(||v| 2 + , ^\ )). When v\ ^ 0, the vector 

{pV 2 ,pV2V,pV2{-\-V\ + _ — )) 

can be expressed as a vector function F(U) of U. Then the system can be written into the form 
in (jl.ip with (t,x) = (x\,X2). By direct calculation, we find that the eigenvalues of this system 
are 

x V2 , , . viv 2 ± ca/|v| 2 - c 2 

Ao = — (repeated), A± = ^ — 5 , 

v\ vf — c z 

and the corresponding eigenvectors are roi,ro2, and r±. Then 

r± • VA± 7^ for all |v| > c, 

r 0i -VA = fori = 1,2. 

Thus, the two characteristic fields corresponding to A± are genuinely nonlinear, while the other 
two fields corresponding to Ao are linearly degenerate. 

An asymptotic geometric optics expansion is of the following form: 

U £ {t,x) = U° + eV £ {t,x), 

where U° is a constant background state in (jl.ip . By a formal derivation of the expansion of 
weakly nonlinear geometric optics for conservation laws (c/. DiPerna-Majda [9]), the expansion 
is expected to be 

n 

U(t, x) = U° + eY, <? U) (e*, x - A^t/^^r^C/ ) + o(e), 
3=1 

and the functions tr^) (r, y) satisfy a decoupled system of scalar conservation laws: 

dr<r® + tfdyia^) 2 = 0, 1 < j < n, 

with the coefficients 

b°j = ■ V 2 / F(C/ )(r,([/°),r J ([/ )) = ^(£7°) ■ VuX(U°). 

Then the genuinely nonlinear condition from fjl -4j) implies that = ^ 7^ which yields that 
equation (|1.5p is the inviscid Burgers equation: 

dra U) + U y (a^f = 0, (1.5) 

while the linearly degenerate condition implies that b® = which yields a linear equation 

d T a^ = 0. (1.6) 

We define 

n 

U £ w {t,x) = U Q + e^<T®(et,x-\ j (V°)t)T j {U Q ). (1.7) 
3=1 
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Then the much simpler function can be used to study the behavior of general entropy 
solutions of nonlinear hyperbolic system of conservation laws, provided that the convergence 
can be rigorously verified. 

An approach has been introduced in Chen-Christoforou-Zhang [21 [3], based on the results 
presented in Bressan [T], to compare the solutions of two different systems, which requires that 
one of them is the standard Riemann semigroup (SRS) while the other is only a global entropy 
solution with bounded variation obtained by the front tracking method. Let T> C L 1 (lR;lR n ) be 
a closed domain. A map S : V x [0, oo[i — > V is a SRS generated by system if the following 
three conditions hold (cf. [I]): 

• Semigroup property: For every U € T> and t, s > 0, 

S U = U, S t S s U = S t+s U; (1.8) 

• Lipschitz continuity: There exist constants L\ and L 2 such that, for all U, V G T> and 

s,t > 0, 

\\S t U - S a V\\ L i < Li||C7 - V\\ L i + L 2 \t - s\; (1.9) 

• Consistency with the Riemann solver: For any piecewise constant initial data U € 
V, there exists 5 > such that, for all t € [0, 5], the trajectory U(t, •) = StU(-) coincides 
with the solution of the Cauchy problem (jl.ip obtained by piecing together the standard 
solutions for the Riemann problems determined by the jumps of U. 

Following [2 [3], in this paper, the general entropy solution U(t,x) of the Cauchy problem 
(jl.ip under consideration is the SRS, which can be constructed by the front tracking method 
(cf. [13]), and the entropy solution for the corresponding asymptotic scalar equation (|1.5p is 
constructed by polygonal approximations, first introduced in Dafermos [7], with initial data: 

a^(0,y) = \ j (U°)-U 1 (y), l<j<n. (1.10) 

We establish the L 1 -estimate between U(t,x) and U^,(t,x) by using both the properties of the 
wave- front tracking algorithm and the standard error formula (cf. [1]): 

PrW( ) - W(T) ht < L [ T ^ S » W ^-Y (T + h)hl <iT, (1.11) 

where L is the Lipschitz constant of the semigroup St, and W(t) is any Lipschitz continuous 
map defined on [0, T]. One of our objectives here is to develop a new approach to provide a 
rigorous mathematical proof of the following theorem. 

Theorem 1.1. (Main Theorem). Let F(U) € C 2 (M n ; W 1 ), and let U 1 ^) be an arbitrary 
function of bounded variation with compact support. Assume that each eigenvalue of the hyper- 
bolic system in (jl.ip has constant multiplicity and its corresponding characteristic field is either 
genuinely nonlinear or linearly degenerate. Consider an entropy solution U £ (t,x) of the Cauchy 
problem (jl.ip . which is a SRS, and the weakly nonlinear geometric optics expansion function 
U^(t,x) defined by (jl.5p ~ (jl.7p and (jl.lOp . Then there exists eq > such that, for all t > and 
e e (0,£ ], 

\\U e {t,-)-U £ w {t,-)\\ L ,<Ce 2 , (1.12) 
for some C > independent of e and t. 

We remark here that this result allows the eigenvalues of the n x n hyperbolic system in 
(jl.ip to have constant multiplicity and the corresponding characteristic fields to be linearly 
degenerate, which answers the open problem posed by Majda in [T2]. In particular, for the nxn 
system in (jl.ip . we obtain that the L 1 — estimate between the entropy solution and the geometric 
optics expansion function is bounded by 0(e 2 ) that is independent of t S [0,oo). 

The proof of Theorem 11.11 is based on our following observation: For a genuinely nonlinear 
system with initial data of compact support, the waves of different families in the solution will be 
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separated each other. This enables us to follow Majda-Rosaales |15] and use the front tracking 
method. We introduce an auxiliary approximate function by adding higher order term of 
to {LTD as 

n 2 n 

V* = U°+e^2a^(et,x- X^t^U ) + j ]T (a®(et,x- X^t)) 2 ^ ) ■ V)r 3 (U°) 
i=i 3 

(1.13) 

and its corresponding more accurate auxiliary approximate function: 

= U° + eJ24 j \zt,x- XjiU^rjiU ) 

j 

+f E WW - ^(U°)t)) 2 (r,(U°) • V)r,([/°) 
+e^E^(x-X j (U°)(t-T )), (1.14) 

where j £ N means that the corresponding j— th characteristic field is genuinely nonlinear, 
while ji € L means that the corresponding characteristic field is linearly degenerate and all {ji} 

together constitute the j-th characteristic field; Furthermore, a)j is given in §2.21 and E{? is 
defined in §5. The novelty here is that the new correction terms are introduced to deal with 
the contact discontinuities. With this key observation, then our approach is to prove the L 1 - 
distance between this auxiliary function and the general entropy solution to system (jl.ip with 
the same initial data is 0(l)e 2 , and finally to employ the L 1 -stability of solutions with respect 
to initial data to establish Theorem II. li The complete proof of Theorem 11.11 will be given in 
§H and ^3 This provides an alternative approach to deal with nonlinear geometric optics for 
hyperbolic systems of conservation laws in (|1.1|) . 

As an example of further applications of this approach, we extend the result to the case when 
the initial data has non-compact support. 

Theorem 1.2. Let F(U) € C 2 (R n ;M n ) andU 1 ^) € BV(R; R n ) n L 1 (R; M. n ) . Assume that each 
real eigenvalue of the hyperbolic system in (|1.1|) has constant multiplicity and its corresponding 
characteristic field is either genuinely nonlinear or linearly degenerate. Consider an entropy 
solution U £ (t,x) of the Cauchy problem (jl.ip . which is the SRS, and the weakly nonlinear geo- 
metric optics expansion U^(t,x) defined by (|1.7p . Then 

sup \\U £ (t, •) - U^(t,-)\\ L i = o(e) whene^O. (1.15) 

0<t<T /e 

For related earlier results in this direction, we refer the reader to DiPerna-Majda [S] for an 
order of 0(et 2 ) for the case of periodic initial data and the same estimate (|1.12|) for initial 
data with compact support for 2x2 genuinely nonlinear and strictly hyperbolic systems of 
conservation laws. For general strictly hyperbolic systems with some kind of periodic properties 
of the initial data for which the resonance phenomena occur, Schochet [16] proved the L 1 - 
estimate of order o(e 2 t); and Cheverry [6] dealt with more general initial data and proved that, 
for all t > 0, 

||C/ e (t, •) — m e (t, •)||l 1 (x) = °{ £ ) f° r an y fixed compact set K (<= K, 

where m e is the corresponding geometric optics expansion. Both of their results allow the 
characteristic fields to be linearly degenerate, but require the hyperbolic system in (jl.ip to be 
strictly hyperbolic. We also refer the reader to Chen-Junca-Rascke [1], Cheverry [5], Gues [10j . 
Hunter-Majda-Rosales [11], Joly-Metivier-Rauch [12], Majda-Rosales [15], and the references 
cited therein for related results. For classical results on the front tracking method and hyperbolic 
systems of conservation laws, see Bressan [1] and Dafermos [8] . 
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2. Front Tracking Schemes and Standard Riemann Semigroups 

In this section, we analyze entropy solutions of hyperbolic systems of conservation laws and 
front tracking algorithms for scalar equations for subsequent development. 

2.1. Existence and Stability of Entropy Solutions. Consider the Riemann problem of 
system (jl.ip with the following initial data: 

f U+, x > 0, 
U(0,x) = { (2.1) 
U-, x < 0, 



where U± are constant vectors. Based on the results in |13j . we have 

Lemma 2.1. Assume that F(U) satisfies the same assumptions as in Theorem 11.11 Then, for 
every compact set K <s there exists 5 > such that, whenever U- € K,\U + — U-\ < 8, 
the Riemann problem above has a unique entropy solution, which consists of n — p + 2 constant 
states, denoted by £/» (i = 0, 1, • • • ,m,m + p, ■ ■ ■ ,n), and n — p + 1 elementary waves (shock 
or rarefaction waves corresponding to the genuinely nonlinear characteristic fields, or contact 
discontinuities to the linear degenerate fields). Moreover, there exists a unique small parameter 
vector • • • , /3 n ) such that 

U = U^, U n = U + = ¥(E7_; ,/3 n ), 

(2.2 j 

U i = ip i (U i - 1 ;Pi), i = 0, 1, - - - ,m,m+p+l,--- ,n, 

and 

U m +p — ip(U m ] /3 m+ i,--- ,/3 m+p ), (2-3) 

where ifr^ij), and are smooth functions with respect to the respective parameter vectors and 
satisfy 

il>i(U-;0) = U-, j^(U-;0) = r, (C/_), i = 0, 1, • • • , m, m + p + 1, • • • , n, 

^-eKer(X(U)I-VF(U)), j = m + 1, • • • , m + p, (2.4) 

V(^-;0,--- ,0) = l7_, U([7_;0,--- ,0)=r,(LT_), j = m + 1, • • • , m + p. 

From Lemma 12.11 one can follow the approach in [I] and [T7] to obtain the existence and 
stability of the unique entropy solution that is the standard Riemann semigroup (SRS). 

Lemma 2.2 (Existence and Stability of SRS). Assume that F{U) satisfies the same assumptions 
as in Theorem W.W Then there is a suitably small 5q > such that, given any U € L 1 (M; W 1 ) with 
TV(U) < 5q, there exists an entropy solution U(t,x) = St(x) by the wave-front tracking method 
or the Glimm scheme. The map S : [0, oo[xV i — > T> satisfies that, for all U ,V G T>, s,t > 0, 

S U = U, S t S s U = S t+s U, (2.5) 
\\StU — S 8 V\\li < L\ \\U — V^l 1 + L2 \t — s\ for some constants L\ and L2, (2.6) 

such that the solution is a SRS. Furthermore, for any other SRS S : T> x [0,oo[i — > D, defined 
on a domain D D T>, we have 

S t U = S t U for any U €V, t> 0. 

In Lemma 12.21 the domain T> is defined by the Glimm functional; see the definition in [T] 
(Chapter 7, pp. 151). 
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2.2. The Wave-Front Tracking Scheme for Scaler Equations. We adopt the following 
scheme and related results from pQ; also see [7j for the details. For our problem, we assume 
that the scalar flux function / is a convex function. For fixed integer v > 1, we consider the 
piecewise constant initial data u taking values within the discrete set 2~ V 7L = {2~ u j : j integer} 
and define f v to be the piecewise affine function that coincides with / at all nodes 2~ v j with j 
integer: 

Ms) = S -^f(2^( 3 + 1)) + 2 ~ U(j ^J ) ~ S f(2^j), s G [2-"j, 2-"{j + 1)]. 

For the scheme, we consider the Cauchy problem: 

d t u + d x f u (u) = (2.7) 

with initial data u. First we consider the Riemann problem with initial data u as in (|2.ip and 
u± E 2~ ur L to obtain the following solutions: 

• Case 1. u_ < u + : We define the increasing sequence of jump speeds as 

, fu(wi) ~ fu(wi-l) , -, 

*i = , l = l,---,q, 

wi - Wl-l 

where {wi}f =1 are the jump points of f v and satisfy wo = U- < w\ < ■ ■ ■ < w q = u + . 
Then we can obtain an entropy solution of the above Riemann problem ()2.7p and (|2.ip 
as follows: 

{U— if x < Ait, 
wi if Xit < x < Aj+it, 1 < I < q - 1, (2.8) 
U + if X > X q t. 

• Case 2. u_ > u+: In this case, we define the shock speed as 

= f v {u + ) - fy{u^) 

u + — n_ 

Thus, we can also obtain an entropy solution as the previous case: 

(4- \ f u ~ i£x < At, , . 

wvt.x) := < .j. . , 2.9) 

v ' ; | u + if x > At. y ' 

Next, consider a more general Cauchy problem for (|2.7p with piecewise constant initial data u, 
taking values within the set 2~^Z. We can construct the solution by solving the corresponding 
Riemann problems so that the total number of interactions is finite and the solution can be 
prolonged for all t > 0. For these solutions u u (t,x), we have the following properties: 

TV{u u (t, •)) < TV(u), \\u v (t, -)||l°° < ||«IU» for any t > 0, (2.10) 

IK(V) -u v {t',-)\\ L i <LTV(u)\t-t'\ for any t, t' > 0, (2.11) 

where L is the Lipschitz constant such that 

\f(w) - f(w')\ < L\w - w'\ for any w,w' G [— M, M]. 

Then, using Helly's theorem, we obtain an entropy solution u = u(t, x) defined for all t > by 
compactness, with 

TV(u(t, •)) < TV(u), \\u(t, -JUioo < ||u|| L oc for any t > 0. (2.12) 

We remark here that, since the scalar flux function in our case in §3- §6 is quadratic, the 
front-tracking algorithm could be applied directly. However, we adopt the piecewise-linear ap- 
proximation in our analysis in §3— §6 so that it is more convenient to compare the solutions to 
present our approach. 
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In this section we compare the Riemann solvers to system (jl.ip and the geometric optics 
expansion defined by (|1.7p with the same initial data. From now on, we denote A^ 1 = Xj(U°) 
and r;- = Tj(U°) through the paper. 

3.1. Comparison for the Genuinely Nonlinear Case. First, we consider some properties 
of approximate solutions to the Burgers equation. 

Lemma 3.1. For some 1 < k < n, assume that, for fixed v, a^ir, y) is an approximate solution 
to the inviscid Burgers equation: 

d T w + -d y (w 2 ) = 0, 

constructed by the front tracking method. At the jump point (tq^q), suppose that cr£_ = 

o-l k \T ,y -), cr^l = o-l k \T ,y +), a^l = af} +a, and (t ,x ) = {^,y + A£3>) is the cor- 
responding point in the (t, x)- coordinates to the point (ro,yo)- Then 

(i) If a < 0, the slope of the discontinuity line of solutions to the inviscid Burgers equation 
in the (t,x) -coordinates is 

S B{ttX) (a) = X° k + ai k le + ^e; (3.1) 

(ii) If a > 0, then a = (o-„2 — 2~ u , and the slope of the discontinuity line of solutions 
to the inviscid Burgers equation in the (t,x) -coordinates is 

Proof. First, we deduce from ^2.21 that the slope of the discontinuity line of solutions to the 
inviscid Burgers equation in the (r, y)~coordinates is 

y\r) = <£> + 1 



from the Rankine-Hugoniot condition when a < 0, and 

V v ) = °v- + g 

from ()2.8p when a > 0. Using the fact that, if (t, x) is a point on the discontinuity line with the 
corresponding point (r, y) in the (r, y)-coordinates, we have 

(y = x-Xlt, ^ (S BM (a) = ^ a<0, 
I T = ^ ^ 1 A2^,(a) = fEff , a>0. 

This completes the proof. □ 
Next, we consider the Riemann problem to system (jl.ip with the corresponding initial data: 

n 

U\ t=tQ =U + eY,4 3) (zt ,x-\°t )r°. (3.3) 

3=1 

Lemma 3.2. Assume that F{U) satisfies all the assumptions as in Theorem W.W the k-th char- 
acteristic field is genuinely nonlinear, and a;? has a jump point as Lemma 13.11 Let (to,xo) 
be the corresponding point to (ro,yo)> U± = U° + eX^=i a ^± r j> an< ^ ^+ = ^K^-j/^lj'" ,fin), 
where = a^l are constant valued functions near (ro,yo) when j ^ k. Then we have 

/3j(a, e) = aedjk + 0(l)<r(max \a^l\ + a) e 2 , 

3 
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where 5jk is the Kronecker delta. 

Proof. Let 9 = ea. Prom the definition, ?7f = f/f. + 9r k . Then 

$(Ut;Pi,--- ,P n )-Ui = 0r° k . 

If 9 = 0, we deduce from <&({/£.; • • • ,/3 n ) = f/f. that ^ |^=o = 0- We differentiate both sides 
of the above equation with respect to 9 to obtain 

/ j ao an L ki 



J 



8^ 89 



From Lemma EH we know 



Then we find 



<9$ 



8$ 

8=0 ~ d $3 



r k {Ut 



/3=0 



89 



e=o 



8=0 



3=1 



= r k (Ul) - £ • V)r fc (£/1)<#2 + 0(l)e 2 

5=1 

= r fc (C/l) + 0(l)max|<7?2|e. 

j 

Taking the dot product both sides of (j3.4|) with lj (£/£.), we obtain 



90 



djk + O(l) max | <T^i | £• 



Then we have 



$(0, e) = S jk 9 + O(l) max \a { Jl\9e + 0{l)9 l 

j 

= 5j k ae + 0(l)a( max \a^l \ + <r)e 2 . 



(3.4) 



□ 



This completes the proof. 

With Lemmas 13. 1H3.2| we have 

Proposition 3.1. Assume that F(U) satisfies the assumptions as in Theorem 11.11 the k-th 
characteristic field is genuinely nonlinear, is a piecewise constant function as Lemma 13.11 

and has a jump point at (to,£/o)j (to, x o) is the corresponding point, and (to,?/o) is not a jump 
(i) 

point of o)j for j ^ k. Then, for every A > 2maxi£n<jif \k\U) with M being the maximum of 
the solution for (|1.1|) . if h is sufficiently small, we have 

\Sh(KJt ,-)) - U*(t + h,x)\dx < Ca(a + max\a { jl\)he 2 , 



(3.5) 



xo—Xh<x<xo+Xh 



where 



U £ w {t,x) := U° + e^al j) (et,x- Ajt)r°, 

3=1 

and C > is a constant independent of e and t. 
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Proof. First, we consider the case a < 0. Then, from Lemma 13.21 we obtain that f3 k < when e 
sufficiently small. Thus, the /c-th wave is a shock. From Lemma 12. f I and the Rankine-Hugoniot 
condition, we have 

= \{r k ■ V)A fc (£/ fc _i), 

where S((3 k ) is the shock speed of the k-th shock with respect to and U^-i is defined as in 
Lemma [2TT1 Therefore, if U- = v _ and U+ = u+ , we have 

S{p k ) = \ k {U k ^) + \/3 k + 0{l)P 2 

= A*(17_) + \fa + 0(l)/3 2 + 0(l)a(a + max,- \a9}\) e 2 , 

by using Lemmas 12. f I and 13.21 and the following estimate: 

\U k - X -U-\ <J2\Uj -Uj-x\ <C^\pj\ = (7a(max|4-| + a)e 2 . (3.6) 

j<k j<k 3 

Using Lemma 13.21 we have 

S{h) = X k (U^) + lae + 0(l)a(a + m a ^\a { Jl\)e 2 

= \ k (U-) + \ae + 0(l)a(<7 + max \a { Jl\)e 2 , 
2 j 

Then we have 

S(f3 k ) - S B ^ x ){a) 

= A fc (C7_) + \ae + 0(l)a(a + max, \a { Jl\)e 2 - (A° + afle + § e) 
= A fc (C/_) - Ag - ia^e + 0(l)ff(<7 + max,- |4 J 2|)e 2 
= + 0(l)ama Xi (|^2| 2 )e 2 + 0(1)<t(ct + max,- |4 J 2|)e 2 

= O(l) max, |cj^_| £■ 

As shown in Figure [H with estimate (|3.6p . and 

- y+| < \ U i ~ U i+i\ < C Y.\^\= ^( max Wu-\ + (3-7) 

j>k j>k 3 



s{p k ) 



Pk=0 



\k{U k - 



dS((3 k ) 



d/3 k 
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we have 



/ . . \S h (U^ u (t ,-))-U^ u (t + h,x)\dx 

J XQ—\h<x<XQ+\h 

n 

= I . \U{t + h,x)-U°-^2ai j \e(t + h),x-X J {t + h))r° j \dx 

J xq— \h<x<XQ+\h j_q 

min{Sfc,S B ( t]3 .)}ft rx +\h 
+ / 

xo-Xh Jmax{Sfe,S B(t)I )}/i' 

n 

x\U(t + h,x)- U° -J2 a uH^(to + h),x- A°(i + h^r^dx 

3=0 

/•max{Sfe,S : s( t x )}h n 

+ / . |f7(i + M) -*7°-^4%(t + /i),z- A J °(t + /i))r°|dx 

0(l)a((T + max \a { Jl\)\he 2 + 0(l)amax \a [ Jl\he 2 

3 3 

O (l)a( a + max \a { Jl\) he 2 . 

3 




Figure 1. Comparison for the genuinely nonlinear case: shock wave 
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FIGURE 2. Comparison for the genuinely nonlinear case: rarefaction wave 



We now consider the case a > 0. As for the case a < 0, we obtain that (5 k > when e 
sufficiently small. Then the k-th wave is a rarefaction wave. From Lemmas 12. II and 13.21 we have 

A fc (E4-i) = \ k (U-) + 0(l)a{a + m^\a { Jl\)e 2 , 

3 

\k{U k ) = X k (U + ) + 0(l)a(a + m a x\a { Jl\)e 2 , 

3 

h(U k ) = \ k (U k ^) + p k + 0(l)a{a + max\al 3 l\)e 2 

= X k (U-) + -ae + 0(l)a(a + ma^\a ( Jl\)e 2 . 
2 j 

These yield 

A fc (l7 fc ) - A fc (Z7fc_i) = \ae + 0(l)a(a + max \a { jl\))e 2 = 0(l)ae, 

6 3 

= X k (U-) + \ae - Ag - a { u % - § e + 0(l)<r(<7 + max, |cr2l)e 2 
= A fc (C/° + e E"=i a?2rj) - A° - a«e + 0(l)a(<r + max, |«r?2 |)e 2 
= 0(l)e E^fcC'i • V)A2<t2 + 0(l)ff£ + 0(l)ff(a + max, |<#2 |)e 2 
= O(l) (<r + max, |cr£2 |)e, 
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and 

\ k {U k ^) - \% x) {a) 

= \ k (U-) - A° - af_e - f e + O(\)o{o + max,- \a^l\)e 2 

= A fe ([/° + e E] =l <#2r?) - A° - age - f e + 0(l)a(a + max, \o® \)e 2 

= 0(1) (cr + max, |cr^|)e, 
where we have used the above estimates and the properties of the centered rarefaction waves: 

{ U- iff <A fc (E/-_), 

[/(*,*) = < [/+ iff > A fc (I7_), 

[ R k (*)(U-) iff = A fc ( J R fe (a)(C/_)) G [A fc (tf_),A fc (E/- + )] s 

when £/± are connected by the fc-rarefaction curve, and R k (o){U-) is a smooth function for a 
and £/_. 

As shown in Figure El similar to the case a < 0, we have 
. \S h (U^(t ,-))-K^t + h,x)\dx 

xo—Xh<x<xo+Xh 

r. n 

= I , \U(t + h,x) -U° -^2a^(e(t + h),x- A°(i + /i))r°|dx 

J xo—Xh<x<xo+Xh j_Q 
r mm{X k (U k ^),X^ } } ™ 

= I ' \U(to + h,x)-U -Y,° { J ) (z(t + h),x-X j (to + h))r°\dx 

Jxo-Xh j_ 
rxo+Xh n 

+ k \U(t + h,x)-U° -^(jW(e(to + /i),^-A J (to + /i))r J |dx 

7max{Ajfe(f7fc),A^ ( ) tjlB) } j=Q 
rmax{X h (U k ),X^l ,} " 

+ / '* \U(to + h,x)- U° -J2°l j) (<t + h),x- Xp + fc))r?|ds 

= 0(l)cj(cj + max|cj^2|)A/ie 2 + 0(1)ct max |cr^2 |/i£ 2 

= 0(l)CT(cr + max|cj^2|)/ie 2 . 

3 

This completes the proof. □ 

3.2. Comparison for the Linearly Degenerate or Constant Multiplicity Case. In this 
case, we first show that the characteristic fields whose eigenvalue has constant multiplicity, bigger 
than one, must be linearly degenerate. 

Lemma 3.3. Assume that F(U) satisfies the assumptions in Theorem and X m +i(U) = 
••• = X m + p (U) for p > 1. Then the characteristic fields of A m +i (U),-- ■ ,X m + p (U) must be 
linearly degenerate. 

Proof. Assume that X m+ i(U) = ■■■ = X m+p (U) = X(U), and the corresponding eigenvectors 
r m+ i(U), • • • , r m+p (U) are linearly independent. Performing tj • V to both sides of the equation 
VF(U)r k (U) = X k (U)r k (U), we have 

V 2 F(r k (U), rj (U)) + VF(rj(U) ■ V)r k (U) = ( Tj (U) ■ V)X k (U) r k (U) + X k (U){rj(U) ■ V)r k (U), 

(3.8) 

where j, k € {m + !,-•• , m + p} and j ^ k. 
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Taking the dot product on both sides of (|3.8p with lj from the left and using the fact that 
lj ■ r k = 5jk, we have 

lj(U) ■ V 2 F(r k (U),rj(U)) + Xj(U) lj(U) ■ (rj(U) ■ V)r k (U) = X k {U)l 3 {U) ■ ( Tj (U) ■ V)r k (U), 
which implies 

lj(U) ■ V 2 F(T h (U),Tj(U)) = if j + k. 

Next, taking the dot product on both sides of (|3.8|) with 1^ from the left and using the fact that 
lj ■ r k = 6 jk , we have 

h(U) ■ V 2 F(r k (U), rj (U)) + X k (U)l k (U) ■ (r y (C/) • V)r k (U) 
= X k (U)l k (U) ■ (tj(U) ■ V)r k (U) + (v k (U) ■ V)X k (U), 

which implies 

v k {U) ■ VX k (U) = l k (U) ■ V 2 F(r k {U), Tj(V)) = 0. 
This completes the proof. □ 

Now we consider the /cth-characteristic field whose eigenvalue has constant multiplicity for 
m + 1 < k < m + p in (jl.2p . In this case, the corresponding scalar case is quite simple, and the 



equations are reduced to be d T o~^ = 0. Thus, ai k \r,y) is a piecewise function and 



(ft) -t ^ 

<7v+ if y > y , 

(ft) - f s 

<r v i if y < yo- 



Then the slope of the discontinuity line in the (t, x) -coordinates is 

(fc ) , _ x - x _ y + (A°r)/e - y - (A°T )/e _ y - y . , . * 

X B(t, x )W ~ JZ^ ~ T/£ _ To/£ " + A ft " V (3-9) 

Next, we consider the Riemann problem to system (jl.ip with the corresponding initial data 
as before. We obtain 



Lemma 3.4. Assume that F(U) satisfies all the assumptions in Theorem W.W m+1 < k < m+p 
as in (fL2D . al k) has a jump at (ro,yo)> o,nd = + °~ ■ Let (to,xo) be the corresponding 
point to (r , y ), = U° + e YJj=i a v± r j> an -d U e + = $(E7I; ft, ■■■ ,/3 n ). Then we have 

f3j(a, e) = 5j k ae + 0(l)a( max |cr^2| + cr)e 2 . 

j 

Proof. Let 9 = ea. From the definition, = Ut + Or®, that is, 

$(£T;ft,--- ,p n )-U- = 6r° j . (3.10) 

If 9 = 0, we deduce from <&(£/£.; ft, • • • , ft) = U- that f3j\e=o = 0. We differentiate both sides 
of ()3.10p with respect to 6 to obtain jjwr^j- = ? k - From Lemma I27L1 we know J^- = 
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(?/£.) . Thus, we have 



a/3, 



/3=0 



2^ dpi 

j 



e=o 86 



0=0 



n 



3=1 



= r fc (C7I) - £^>° • V)r fc (C/I)a^ + 0(l)e 2 

3=1 

= r fc (C71) + 0(l)max|a?2|e. (3.11) 

3 

Then we take the dot product on both sides of ()3. 1 1 j) with lj(Ut) to obtain 



Ml 
86 



= 5jk + 0(1) max | (7^2 |e, 
6»=o i 



and then 

e) = 5 jk 6 + O(l) max \a [ Jl\6e + O(l)0 2 

3 

= 5 jk ae + 0(l)a( max \a^l \ + crje 2 . 

This completes the proof. □ 

With Lemmas I3.3fl3.4t we have 
Proposition 3.2. Assume that F(U) satisfies the assumptions in Theorem II .11 either the k-th 

(k) 

characteristic field is linearly degenerate or its eigenvalue has constant multiplicity, at, is a 
piecewise constant function which satisfies crS = + a and has a jump point at (tq^q), 

(to,xo) is the corresponding point, and (ro,yo) not a jump point of o~„ f or j 7^ k. Then, 
for every A > 2maxi u i<j\ ! f X k (u) with M being the maximum of the solution to (jl.ip . if h is 
sufficiently small, we have 



I 

Jxn 



\S h (U^ u (t , •)) - U^ v (t + h,-)\dx <Ca(a + max \a { Jl\)he 2 , 

I xo—Xh<x<xo+Xh 3 

where v is defined in Lemma [3. 11 and C > is a constant independent of h and e. 

Proof. The proof of this lemma is similar to the one for Lemma [3. 11 especially for the shock case. 
The only change is about the shock speed S(f3 k ). This can be calculated directly as follows: 

fiW = \ k (U k -i) = \ k (U-) + 0(l)a(a + max \a^_\)e 2 = X° k + 0(1) max \a ( Jl\ e, 



and then 

A. |0 

3 

r(?) I 



^fc)-^)^) = A^ + 0(l)max|4 J 2|e-A° 



0(1) max |<rXl|e. 

3 



With this, we can obtain the desired estimate. □ 
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4. Proof of Theorem [HQ Part I: < t < T 

We now prove Theorem l 1 . 1 1 when t € [0, To], where To is denned as follows: Since supp(cr^ | r — q) 
is a compact set for each j (due to the compactness of the support of U 1 ) and the speed A^ ^ \® 
for k j, there exists To > 0, independent of e € (0, Eq], such that, for t > To, the compact sets 

Kj(t) = {x € R : a j (et,x - Ajt) ^ 0}, j = 1,2, ••• ,n, 

are disjoint, that is, for t >Tq, 

Kj(t)nK k {t) = 4> xij + k. 

In order to establish Theorem 11.11 for < t < To, a technical lemma is stated here to restrict 
the integral interval for the semigroup estimates to be finite. 

Lemma 4.1. Let S : T> x [0, oo) i— >■ T> be the semigroup generated by system (jl.ip . and let A 
be an upper bound for all the wave speeds, where F{U) satisfies all the assumptions in Theorem 
11.11 Given any interval Lq := [a,b], define 

r c , c i b — a 

I t :=[a + Xt,b- At , t< — — . 

2A 

Then, for every Lipschitz continuous map W : [0, T] i— > V, 

\\W(t) - S t W(0)\\ L i (It) <L J jlimmf ^ j dr. (4.1) 

Proof. First, for every U, V G T> and t > 0, we have 

b-Xt rb 

\S t U(x) - S t V(x)\dx < L / \U(x) -V(x)\dx. 

a+\t Ja 

This can be obtained by using the following facts: 

(i) If two initial conditions U, V € T> coincide on (a, 6), then StU(x) = SfV(x) for all x € 
(a + At, 6 -At); 

(ii) HtJeV, then [7 x [a ,6] e T>. 

With this estimate, the other steps for the proof is the same as the ones for ([l.lip in [1] 
(Theorem 2.9). □ 

Now we prove Theorem II. II for < t < To. 

Proof of Theorem 11.11 for < t < To- We divide the proof into three steps. The constant 
C > below is a universal bound independent of the parameters (j, v, e, r, t, J) in the proof. 

Step 1. Approximation by piecewise constant functions. Since U 1 £ .BV, we can construct 
piecewise constant functions a^(0, ?/), which satisfy 

TV(<tP{0,-)) <TV(ad\0,-)), 

\\al j) (0, •) - a<i\0, OHoc + 11^(0, ■) - <t®(0, -)\\^ < e, 
where <r®(0,y) = lj(U°) ■ U l {y). Then 

TV(a®(p, •)) < TV(a U) (0, •)) < CTl^f/ 1 ) 

and 

ll^'Ho, OIU- < 211^(0, <c\\u l {-)\\ L ~. 

In addition, the jump points of each are finite, and the range of those functions is 
contained in the set 2 _I/ Z, for some fixed integer v. For this initial data, we can construct an 
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approximate solution Ov (r, y) by the method introduced in £ )2.21 Then the solution cr£ (t, y) 
satisfies 

TV(a^(r, •)) < CW(<#)(0, •)) < CTV(U l {-)), (4.2) 

lk^(T, -)|U<- < C \\ajp(0, < CH^COIIl-, (4.3) 

\\aW(T,-)-aP(s,-)\\ L i<L\T-s\. (4.4) 

Then, from Helly's theorem, for every r, there exists a subsequence of functions (still denoted 
as) Gv (r, y) so that 

a® (r, y) — ► o-C?) ( y ) for a.e. y G R, 

and 

< CII^OII^, TV(a^{-)) < CTV(U\-)). 

By standard argument, there exists a subsequence cr^v such that cr^v (t, •) — > o"( J )(r, •) pointwise, 
which is also in Lj oc (R;R n ), at every rational time r > 0. Define 



a 



(r,y) = lim a (j) (r m ,y), 



where r m is a rational time and r m -4 r. Then, by continuity of time in (j4.4f) . we conclude that 
<7^ -4 in ^([0,00) x R;R n ), and a® satisfies 

/oo 

|<t^(t, 1/) — a^(s, y)\dy < L\t — s\ for all r, s > 0, 
-00 

3V(ff W (T, •)) < CTV^C/ 1 ), H^Ct, OIIl- < C\\U\')\\ L «, for all r > 0. 

Step 2. Estimate of the term \\ShU^ v (t, ■) — v (t + h, OIIl 1 ^) for h small enough, where 
I:= [a,b], 

n 

3=1 

and (i, x) is not a point of interaction of (et, x — \®t) for all x £l and j = 1, • • • , n. 

First, consider the case that (t, xq) is a jump point of one and only one function of |<7^ (t, 
say o~v\r,y). From we can derive the following estimate: 

ll'S'/i(^S, v (*» 0) - C^ >v (* + ^, OIUiW ^ - ^S|(max 1^21 + max lo^D/^e 2 , (4.5) 

3 3 

where cr£2 := lim x ^,. Xo ± o~l .{et,x — A^t), when h is so small that there is no interaction, and I 
contains only the wave fronts generated at point (t, xq). 

We now consider the general case. Without loss of generality, assume that (t, xq) is a jump 
point of ai^ and . 

Claim. Assume that F(U) satisfies all the assumptions in Theorem II. 1L j and k are two 
distinct integers, (ro, yo) is a jump point of and only of o~„ and al k \ Let 



and 



a { jl= lim a^(T Q ,y) 



a£U lim ai k \r ,y), 
y->yo± 
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(to,xo) is the corresponding point to (ro,yo)> 

n 

i=i 

Then we have 

Pm — [o m j \cr y+ - o- v _ I + d mk \a v+ - a v _ \ ) e 

+0(l)(kS -<&-\ + -offi|)(m«|a2l +max|aS|), 2 . 

This property can be proved by combining the properties in the proof of Lemma [3 .21 and the 
ideas for the proof of Lemma 13.41 We omit the details. 
With this, we can show 

\\S h (U^ u (t,-))-U^(t + h,-)\\ L t (I) 

< C{\*$-*21\ + |ag - a«|)(in«|a£| +ma X \*%\)he\ (4.6) 

j j 

when h is so small that there is no interaction and / contains only the wave fronts generated at 
the point (t, Xq). The proof is similar to the one in Propositions 13. ll and 13.21 the difference here 
is that we should divide the interval into three parts: the additional part is near the jth-wave. 
Thus, we estimate not only — SB(t,x)(°~) or ^k{Uk) — A^i^o") as * n P r °P os iti° n 13.11 or 

E2 but also 5(^0 - S BM (a) or \j(Uj) - \% iX) (a). 

With (|4.5p and ()4.6p . we can derive an estimate for arbitrary finite interval I by dividing 
the interval into some subintervals, every one of which contains only one jump point. Since the 
number of these subsets is finite, we can sum together to obtain 

\\Sh(UZ, v (t,-)) -U^ u (t + h,-)\\ LHI) < CY,TV {t)X)mxI (a^(et,y))m^\a^\he 2 , (4.7) 

where (r, y) = (et, x — A°t). 

Step 3. Completion of the proof. For an arbitrary finite interval J = [a,b], let I = [a — 
Xt,b + At]. Then, from the definition in Lemma 14. 1\ we have It = J- For any initial data 
U l {x) € BV n L 1 , a® and U*(t,x) are defined from (fT3|) . ([LTD . and (fTTTOl) . Then, following 
Step 1, we construct piecewise constant functions a„(T,y) so that 

\\u%t, •) - i%(t, OIUhj) < ||f/ £ (t, •) - t£(t, -)IUi(j) + \\K(t, •) - ^(t, -)IUi(j) 

+ ||[/-(t ) .)-C/^( i ,.)|| il(J) 

= h + h + h. 

For Ii, from Lemma 12.21 



/! = llSt^ + e^OJ-St^ + eJ^^^OrJJIkifj) (4.8) 
< L||[/ + £ C/ 1 (-)-(t/ + ^^)(0,-)r°)|| il (4.9) 



< Cj2h U) &-)-4 j \0,-)\\ Ll . (4.10) 
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Hence, I\ —> when v — > oo. 

For I2, it is reduced to estimate ||oir(r, •) — a^\r, OIIl 1 ^)- Prom Step 1, we can select a 
rational time sequence {r m }^ =1 and a subsequence (still denoted) {u} such that r m — > r when 
m — > 00 and, for all m, (Tv(T m ,y) — > o"W(r m ,?/) for a.e. y € M, when i/ — )■ 00. Then 

Ikv '(Tm, •) - o- (i) (r m , OIU^j) wnen ^ 00. 

Using the Lipschitz continuity (|1.9p . we obtain that J2 — >• when v — >• 00. 
For I3, from Lemma 14. 11 we obtain 

^3 = ll^(i,-)-^(t,-)||LM/ t ) 

. . ll^(^ + ^-)-^.(^(^-))||Li(/ s+fe ) , 

< L I hmmi - — ds 

h 



C £ 2 /* E"=i ( £Sj •)) max, \\(Ti 3 \es, , 

< L hmmf ds 

Iq h^o+ h 



< Ce 2 [ (TV{U r )) 2 ds 
Jo 

= C (TV(U l )) 2 te 2 , 



where we have used the fact that TV(a { J\et, •)) < CTV(a ( J\o, •)) for all t > 0. With these 
three estimates together, passing limit v — > 00, we obtain that, for any arbitrary finite interval 
J := [a, b], 

\\U"(t, ■) - U^t, OlkV) ^ C (TViU^fte 2 . 
Finally, let a — s> —00 and b —> 00, we have 

\\U%t,-) -U^(t,-)\\ L1 <C (TViU^Se 2 . (4.11) 

Therefore, for < t < To, we obtain 

\\U e (t,-) - U^t,-)\\ L i <CT e 2 . (4.12) 

□ 

As a direct corollary, we conclude 

Corollary 4.1. Assume that F satisfies the assumptions in Theorem 11.11 and U 1 £ n 

L 1 (M;IR n ) (not necessarily with compact support). Consider an entropy solution U e (t,x) of 
(jl.ip . which is the SRS, and the weakly nonlinear geometric optics expansion U^(t,x) defined 
by (|1.5p . (jl.7p . and (jl.lOp . Then there exists > siic/i i/iai, u;/ien e € (0,£o], we /mue 

||?7 e (t, •) — U^(t, -)\\li < Cte 2 /or some constant C > independent of e and t. (4.13) 

However, estimate (|4.13p is not really strong enough to justify the weakly nonlinear geometric 
optics approximation as noted in Schochet [16] . In §6, we further develop the approach to 
improve (|4.13p into the stronger error estimate (jl. 15[) in Theorem 1.2, which is strong enough 
indeed to justify the approximation even for the case when U 1 G BV n L X {R\ W 1 ) without 
compact support. 

5. Proof of Theorem 11.11 Part II: t > T . 
Notice from SQjthat, when t > To, the compact sets Kj{t),j = 1, 2, • • • , n, are disjoint: 

K j (t)nK k {t) = <f> Hj^k. 



NONLINEAR GEOMETRIC OPTICS FOR HYPERBOLIC SYSTEMS OF CONSERVATION LAWS 19 

In this section, we give a careful calculation based on this and complete the proof of Theorem 
11.11 The constant C > below is the universal constant, independent of (e, v, h, t, r) in the 
proof. 

We first refine expansion (|1.7p or (|3.3|) by the following corresponding approximation: 
V? = U + e^a^\et,x-X 3 t)r J + £ ^^(a^(et,x-X J t))\r° r V)r J 

+e 2 ^2E^(x-X° 1 (t-T )), (5.1) 

where j G N means that the corresponding j— th characteristic field is genuinely nonlinear, 
while ji G L means that the corresponding characteristic field is linearly degenerate and all {ji} 
together constitute the j-th. characteristic field. Furthermore, 



and Wv is defined as follows. First set [a, b] := supp x a^^ (eTq, x — X®Tq) and let Xj i G [a, b) be 
a point at which cr^(e?o, x — A^To) has a jump. Denote the set of these jump points together 

( 7 ) 

by Jj and then define the piecewise constant function W{? with the same jump points in Jj. 
At each jump point Xj. G Jj, if 

^S. ) -o-?i ) )rg = *(^ 1 ,-.-,/3 B ), (5.2) 

and the multiplicity of the j-th characteristic field is m, then the difference of Wu at this jump 
point is defined by 

W$ - W® := $(0, • • 0, fa , ■ ■ -P jm , 0, ■ ■ ■, 0). (5.3) 

For the point x < a or x > b, we define Wv — U$. Then we have the property that Kj(t) = 

swpY> x {Wu\x — \j(t — To)) — U }. From the proof of Lemma [531 later, we will see that the error 

(i) 

term satisfies that \E{, \ < C. 

Now we consider the case j G N. As the bootstraps in £J2J we have the following lemma. 
Lemma 5.1. Suppose that the assumptions in Lemma \3.2\ hold. Let 

K ± := [7° + ea?M + j ■ V)r° 

with $(VjJ_;/3i, • • -,/3 n ) = V„ + - Then we have 

= aeb k3 + 0(l)e 3 {(a^lfa + a ( u % 2 + a 3 ) , (5.4) 

where a = a^2 — o~^]_, and 5^ is the Kronecker delta. 
Proof. Let 9 = ea. Then 

V u % - V°_ = Ovl + ^(r° • V)r° + eOa^l ■ V)r°. 
Thus we have the following equation: 

$(Kf_; /3 U ■ ■ n ) = + 9vl + + eOa^l) (r° • V)r°. (5.5) 

Clearly, by the implicit function theorem, we obtain a unique solution (/3i, • • • , /5 n )(#) of (|5.5p 
such that f3j(0) = 0, j = 1, • • • , n. 
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Differentiating (|5.5j) with respect to 9 and then letting 9 = 0, we have 

E ^-^■/•■■ ft -> g| U _ rj + • V)rJ + . V)r?| fe , (5.6) 

Since Ki D Kj = (j) Hi j , we find that uV- 7 = if i 7^ j. Using this fact to calculate (|5.6p . we 
have 

E Sl„ - E^-)§U = r 4 ra + oax^jv 

which yields 



a& 9e ^ JV "- y a# 



^| e=0 = ^ + O(l)(e ) )V. (5.7) 

Next, taking twice derivatives on both sides of (j 5 . 5 j) with respect to and then letting 9 = 0, 
then 



4^ a^a^- a# a# le=0 ^ a^- a# 2 ]e=0 K k J k ' 

% ->3 J 



Thus, we have 



Noting that (r fc (V^_) • V)r fc (T/ J f_) = (rg • V)r fc (C/°) + 0(l)tr2e, we obtain 



a# 2 



9=0 = 0(l)a™e. (5.8) 



Combining (|5.7p - ([5.8p with the Taylor expansion, we finally have 

m =m + ^\e=o d + ^-\e=o d2 + °w 

=ae5 kj + 0(l){(a k u _) 2 a + af_a 2 + a 3 ) e 3 . 

This completes the proof. □ 

With this better estimate, as in £j3l we further have 

Lemma 5.2. Suppose that the assumptions in Theorem W.W hold. Then there exists ho > such 
that, when h G (0, ho], for t > Tq, we have 

\S h (VJ(t ,x))-V„ £ (t + h,x)\dx 

xo—\h<x<xo+\h 

< C7(max|^_|) 2 |a|/ie 3 + C2~ u \a\he 2 . 

Proof. We divide the proof into two cases. 

Case 1 (shock case): a < 0. In this case, by Lemma loTTl f3 k < 0. Denote S k the k-shock 
speed of j3 k in the (t, x)-coordinates and V k = $(VJ_; fii, ■ ■ -,f3 k ,0, • • -,0). Then 

S k = h(Vk-i) + \h + 0(l)\f3 k \ 2 = X° k + al% + l -ae + 0(l)(a 2 + e\a\) e 2 . 

On the other hand, the shock speed of the Burgers equation of a k in the (t, x)-coordinates is 

S Bk = X° k + 4% + l -ae. 

Then 

S k -S Bk = 0(l)(a 2 + \a\\a (k ]\)e 2 . 
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Meanwhile, with the fact that \Vk-i — V n \ + \V k — Vq\ = 0(l)\a\e, we deduce the core estimate 
for the shock case: 

\S k - S Bk \{\V k -t - VU\ + \V k - V U E _\) = 0(l)(\a\ 3 + \a\ 2 \4 k 2\)e 3 . 

Then, with this estimate in hand, following the bootstrap in the proof of Proposition 13.11 step 
by step, using the results in Lemma 15. 11 and performing cumbersome calculations, we have 



f . \S h (V u e (t,x)) -V u £ (t + h,x)\dx <C(m<xx\a*_\) 2 \a\he 3 . (5.9) 

J XQ—\h<x<xo J t-\h 



Case 2 (rarefaction case): a > 0. From the scheme for the scalar equation in H2.21 or by 
Lemma [XU we conclude that the strength of discontinuity satisfies the following property: 

a < C2~ v . 

First, from Lemma 13. 1\ the speed of rarefaction front of the Burgers equation in the (t, x)— 
coordinates is 

On the other hand, the speed of characteristics is 



A fc (F fc _!) = Afc(K-) + 0(1) + ai k 2a 2 + a 3 )e 3 , 

Xk(V k ) = X k (V u \) + 0(l){(ai k 2) 2 a + a^la 2 + a 3 )e 3 , 
\k{V k ) - \ k {V k -x) = 0(l)ae + 0{l){(a {k 2) 2 <J + a { u % 2 + a 3 )e 3 . 
Thus, by following exactly the proof of Proposition 13.11 direct computation yields that 

\S h (V°(t ,x)) - V^(t + h,x)\dx < Ca 2 he 2 < C2~ u \a\he 2 . (5.10) 



/ 



' xo—Xh<x<xo+Xh 

Combining (|5.9p with ()5.10|) . we arrive at the result. □ 

Next, consider the linearly degenerate case, namely, {ii}^ C L. By construction, all the 
values of Wu he on the curve Sj(Uo), i.e. the integral surface of the vector-fields spanned by 
rj i passing through Uq, while A° is the speed. Here W„ is introduced below (|5.ip to replace 
=1 Ov^T^j.- Thus, Wy is exactly a Riemann solution. Then we have 

Lemma 5.3. Suppose that the assumptions in Theorem W.W hold. Then there exists ho > such 
that, when h € (0, ho], for t > Tq, we have 

[ _ _ \S h (V u e (t,x))-V u E (t + h,x)\dx = 0. 

J xo—Xh<x<xo+^h 

With Lemmas 15.21 and 15.31 in hand, summing up the above estimates and using (jl.lip . we 
have 

Lemma 5.4. Suppose that the assumptions in Theorem W.W hold. Then 

\\S(t - T )V £ (T ) - V £ (t)\\ L1(R) < Ce 2 when t > T . 

Proof. We first use the standard error formula (jl.lip and then let h small enough such that 
S(h)VJ is the solution obtained by piecing together the standard entropy solutions of the Rie- 
mann problems determined by the jumps of V„. Then we can use Lemmas 15.21 and 15.31 to 
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obtain 



\\s(t-T )v;(T )-v;(t)\\ L i 



J Tq h^o+ h 



where is the sum of all the jumps at t = z + h with jump strength a. Then, passing the 
limit v — > oo, we have 

\\S(t-T )V E (T )-V%t)\\ L i<C ^ + £ 3/2^^ < Ce 2 5 

where we have used the fact that, for genuinely nonlinear scalar conservation laws, if the initial 
data no has compact support and satisfies ||i*o ||oo — M, then the solution u(t,x) satisfies 

TV(u(t,-)) < Ct - *, 

where the constant C depends only on uq. This completes the proof. □ 



Since the estimates for the auxiliary function have been established, we are now at the stage 
to consider the estimates for the geometric optic expansion function U^. 



Lemma 5.5. Suppose that the assumptions in Theorem W.W hold. Then 

\\S(t - T )UUT ) - K(t)\\ L i {R) <Ce 2 when t > T . 
Proof. This is in fact a simple corollary of Lemma 15.41 Notice that 

K, v (t) =V°{t) - \e 2 £ {*®(et,x - A°t)) 2 (r0 • V)r° 

jeN 

-e 2 ]T4%-A°(i-To)). 

From Lemma l3~4"l and the fact that the discontinuities of Wv and U° + e Y^i=i a » T \ snare the 
same speed A°, we notice that, for any x E Kj(To), x — \*j(t — To) G Kj(t) and 

afreet, x - \°jT - X°j(t - T )) = a^(sT ,x - A°T ), 
wP(x-X° j (t-T )) = wP(x). 
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Thus, we have 

|£0)(*-A°(i-To))| 

1 m 

= -\U° + eY,4 M) (eT ,x- A°T )r° - W®(x)\ 



i=i 



(7) 

where crj € Jj means that there is a jump point of (J{, at which the strength is <7j. 



Thus, by Lemma [ 
||S(i-ro)^(r )-^(i)||Li 

< ||5(t - r )[/^(r ) - S(t - T )V*(T )\\ L i + \\S(t - T )V?(T ) - V?(t)\\iA 

+ \K £ (t)-u^ u (t)\\ L1 

< L\\V U £ (T ) - U^ V (T )\\ LHR) + \\S(t - T )K £ (T ) - ^f(t)|| £ i + \\V u £ (t) - U^(t)\\ Ll 

< \\S(t - T )V*(T ) - V*(t)\\ Ll +Ce 2 Y, 114%^, • - \%)\\ L » ll4%?o, ■ - AjTo)^ 

+ Ce 2 J2 \U j \et, ■ - A^)|| LOO • - A?t)llii + Ce 2 £ |sup Pa! <#>| 

< \\S(t - T )V U £ (T ) - V u £ (t)\\ L i + Ce 2 . 

Then, passing the limit v — > 00 and using Lemma 15.4^ we have 

\\S(t-T )UUT )-K(t)\\ Ll <Ce 2 , 

which completes the proof. 

With Lemmas I5.1H5.5( we can now complete the proof of Theorem 11.11 
Proof of Theorem 11.11 Notice that 

\\s(t)uuo)-K(t)\\ Ll 

< \\S(t)U*(0) - S(t - T )U^T )\\ L1 + \\S(t - T )U^(T ) - U^(t)\\iA 

= h + h- 
Using (|4.12p in we have 

h < \\S(t - T )S(T )U^(0) - S(t - T )U°(T )\\ L1 

<L\\S(T )UUO)-K(T )\\ L1 
< CT e 2 . 

By Lemma |5.5| 

h < Ce 2 . 
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Therefore, we conclude 

\\U s (t,-)-U^(t,-)\\ L i <Ce 2 . 
This completes the proof. □ 

6. Proof of Theorem 11.21 

Following the approach developed in §3-§5, we can extend the result even for BVdL 1 initial 
data, Theorem 11.21 

We first define the error terms E v (t, x) for each time t. For any jump point x of o~l k \st, x — 
X° k t), the jump of E v at x is 

E u (t,x+)-E u (t,x-) = £>£V,z -^)(^V)r° 

+^(ri5) 2 -(^) 2 )(r° fc -VK. (6.1) 

From the properties of a^{et,x - A°t), we know that E„(t,-) € L 1 (R; R n ) n BV(R;R n ) 
and, near the jump point (t,x), the speed of discontinuity of E u (t,x) is the same as the one of 
o~i k \et,x — A^t). In order to use the standard Riemann semigroup, we need to show that the 
error terms E u (t, •) are Lipschitz-continuous with respect to t. Indeed, we have 

\\E u (t,-)-E u (s,-)\\ L i 

= 0(l)||o"^ ||oo(t — s) • [maximum speed] • [total strength of all wave-fronts of al(es, •)] 
< C(t - s) 

for some constant C independent on v. Then, as before, we modify expansion (|1.7|) by the 
following corresponding approximation: 

= U° + e^a { J\et,x- A^)r° + E u (t,x)e 2 . (6.2) 
i 

Then, as the bootstraps in Sj5l we first have the following lemma. 

Lemma 6.1. Suppose that the assumptions in Theorem \1.2\ hold and, at point x, a^ k \et,x — \t) 
has a jump. Let 

Vf± := U° + + E v± (t, x)e 2 

with $(VJ_;Pi,- ■ -,/?„) = V u £ + . Then 

= asS kj + 0(l)e 3 £ ((<x« fa + a ( j ] _a 2 + a 3 ), (6.3) 

i 

where a = a„2 — cr^}, and 5 k j is the Kronecker delta. 



Proof. Let 9 = eg. Then 

V u % - V U L = 9r° k + ^(r° ■ V)r° + e9 £ a£(r° • V)r°. 

i 



Thus, we have the following equation: 

32 



<&(K e -;/3i, • • ■ 1 M = V^+9vl + -{vl-V)vl + e9Y J ^ ) M- v yk- (6-4) 

3 



Clearly, by the implicit function theorem, there exists a unique solution (j3i, ■ ■ ■ ,(3 n )(9) of (fBT 
such that (3j(0) = 0, j = 1, ■ ■ ■ , n. 
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Differentiating ()6.4|) with respect to 6 and letting 6 = 0, we have 

9$(V*_; ft, • • ■, ft) 9ft 



E 



0=0 



00j 89 
r»+^ OT «(r°.V)r»+0(r».V)r» 



0=0' 



(6.5) 



Notice that r° k + e Y^j ^- ( r °j ' V K = r i (K- ) + 0(1) e J2j 1 • Then we obtain from JB3]) that 



Cj)i 



E 



a$ a/3j 



a^ a# 



9/3. 



E r ,(^)^Lo = r fc (^) + o(i)^ 2 E(^-) 2 ' 



0=0 



00 



which yields 



86 



__ = 8 kj + O(l)e 2 J2(^) 2 



Next, takeing the twice derivatives on (|6.4|) with respect to and letting 6 = yield 



E 



d 2 $ 0ft a/3,- 



a/Wj a0 a0 



+E 



a$ a 2 /3 



Thus, we have 



e=o ^ a/3,- a0 2 



(r fc • V)r fe (^_) + 0(l)e 2 ^(a2) 2 + £ r,^ 



0=0 



(r* • V)r°. 



a 2 /?,- 



a0 2 



0=0 



(**-V)rg. 



Noticing that (^(V^) • V)r fc = (rg • V)r£ + 0(l)ecr^l, we have 

a 2 /3,- 



ae 2 



0=0 



(i) 



With all of these formulas above, by the Taylor expansion, we finally have 

0& , „ , a 2 /3,- 



=/3 y (0) + -^| e=o + ^|, =0 ^ + 0(1)0" 



□ 



(fe) 



It completes the proof of this lemma. 

With this better estimate, as in §5\ we have 

Lemma 6.2. Suppose that the assumptions in Theorem \1.2\ hold, and the jump strength of a 
is o~ k . Then there exists ho > such that, when < h < ho, we have 

\S h (V*(to,x))-V°(to + h,x)\dx 

XQ—\h<x<xa+\h 

<C £ \a { Jl\\a k \he 2 +C(^\a { jl\) 2 \a k \he 3 + C2^\a k \he 2 . 
A^A° j 

Proof. The proof is divided into three cases: the contact discontinuity case, rarefaction case, 
and shock case. 

For the rarefaction case, the proof is the same as the one in Sj5l which needs use the crucial 
property of the wave strength in the scheme that a < C2~ u . Thus, it suffices to consider the 
other two cases. 
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Consider a k < for the shock. By Lemma 16.11 j3 k < 0. Denote Sk the k-shock speed of /3k 
in the (t, x)-coordinates and V k = Q(VJ_; f3\, ■ ■ -,f3 k ,0, • • -,0). Then 

S k = A fc (T4_ 1 ) + i/5 fc + 0(l)|/3 fc | 2 



Afc + ea^l + a -f + 0(l)e £ af_ + 0(l)e 2 (a| + e\a k \). 



3+k 

On the other hand, the shock speed of the Burgers equation of o~v in the (i, ^-coordinates is 

n (k) 1 
SBk = Afc + cr^e + -<j fc £. 

Then 

S fc - 5 Bfc = 0{l)sJ2 + + I^I^De 2 . 

Meanwhile, with the fact that |Vfc_i — V n \ + |Vfe — Vol = 0(l)|<Jfc|e, we deduce the core estimate 
for this case: 

\S k - S Bk \{\V k -i - KVI + \V k - V; e _|) = 0(l)e 2 + 0(l)(\a k \ 3 + \a k \ 2 \ai k 2\)e 3 . 

j+k 

Consider the case of contact discontinuity whose corresponding characteristic field has con- 
stant multiplicity: Vj • VX k = if Xj = X k . Then the kth-shock speed S k of j3 k in the (t, x)- 
coordinates is 

S k = X k (V k _ 1 ) = X° k + 0(l)e Y, °u-- 

On the other hand, the shock speed of the Burgers equation of a k in the (i, x)-coordinates is 

Sbu = X k . 

Then we can obtain the same core estimate for this case as for the shock case. With these 
estimates in hand, following the bootstrap in £JS]step by step, we have 

f _ _ \S h (V u £ (t ,x))-V u £ (t + h,x)\dx 

<C Y Wi 3 l\\a k \he 2 + C (YWi 3 l\) 2 Wk\he 3 . 

Finally, combining this inequality with the inequality from the rarefaction wave case, we 
completes the proof. □ 

Then, summing up the above estimates, we have 
Lemma 6.3. Suppose that the assumptions in Theorem \1.2\ hold. Then, for t > 0, we have 
\\S(t)V £ (0)-V £ (t)\\ L1{R) 

<C(et + 2^t)e 2 + Ce 2 £ £ \a^(es,x - X jS )\\a k (x)\)ds, 

J ° k xeMa)\0jt\0 

where x G J k (s) means that, at the point (es,x — X k s), a„ has a jump with the corresponding 
jump strength o~ k (x). 
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This can be achieved by first using the standard semigroup error formula ([l.lip . then letting 
h small enough such that S(h)VJ is the solution obtained by piecing together the standard 
entropy solutions of the Riemann problems determined by the jumps of VJ, and finally using 
Lemma 16.21 

Since the estimates for the auxiliary function have been established, we are now at the stage 
to consider the estimates for the geometric optic expansion function U^. 

Lemma 6.4. Suppose that the assumptions in Theorem 11.21 hold. Let I m = (m,m + 2). Then, 
for t > 0, we have 

\\s(t)um-uut)\\LHK) 

< C(e 2 + eH) + Ce 2 f £ £ \\a^\\ Loa{Im) TV xeIm (a^)ds. 
Proof. This is a simple corollary of Lemma 16.31 Notice that 

and that E„(t, •) G W(M;M n ) n L 1 (R;R n ). Thus, by Lemma we have 

\\s(t)K iU (o)-u^(t)\\ L i {R) 

< \\S(t)U^(0) - S(t)VZ(0)\\ LHR) + \\S(t)Vj(0) - Kf(t)|Ui (R) 

+\K(t)-K >1/ (t)\\ LHR) 

< L\\V u %0) - C/^(0)|| x i (R) + C (st + 2-H)e 2 + \\Vj(t) - K >v (t)\\ Ll{R) 
+Ce2 /'El E E \^\es,x-X jS )\\a k (x)\)ds 

<C(l + et + 2-"(t-T ))e 2 + Ce 2 fj^i E ^ Wi j) (es,x - X jS )\\a k (x)\)ds. 

Jo k xej k {8)\°?\° 

Then, passing the limit v — > oo, we have 

||5(t)C^(0)-C^(t)|| Ll(R) <C(e 2 + eH) + Ce 2 E lk (i) lko 0( 7 m) T14 67m (a«)d S , 

which completes the proof. □ 

With Lemmas I6.1H6.41 we now prove Theorem 11.21 
Proof of Theorem 11.21 It suffices to prove 

lim sup ef £ £ \\a^\\ Loo{Im) TV xeIm {a^)ds^0, 

£ ^°0<t<T /e JO m&x o, x o 

j'k 
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where I m = (m, m + 2). Let r = es. Then 

lim sup e [ E \W ij) \\L^(i m )TV x€lm {a^)d S 



0<t<T / 

-T /e 



E E k {j) h<~ { i m) TV xeIm (<r {k) )te 

/■To A Q >0 

= SS/ E E lk°' ) (r,--fr)||^ (7m) rT4. 67m ( C TW(T,.-^r))dr. 
- 70 meZA^A° 

Since a 3 (r,y) G .BV n L 1 and A° ^ A°, then, for any given r G (0, To) and any 1 < j < n, we 
have 

lim o- (j) (r,y) = 0, lim TV,^ y)u i yoo )(a^) = 0. 

Then, for any ei > 0, there exists a constant Mj > such that, for any \y\ > Mj, \a^(r, y)\ < e\. 
Next, let e > small enough such that, for any A° ^ A°, 

TV" A q_ A o A q_ A o (cr(r, •)) < ei- 

(_ Mj -2+ -L— *r - +2+ -J--JL T ) 

Thus, for any e > small enough, set 

X°- \°- 

J(j, r, e) := {m G Z | I m C (-oo, -Mj + -^-r) U (Mj + -^-r, oo)}. 



Then we have 



E E ll- ' ) (^--^r)|| ioo(/m) TK e , m (^)(r,.-^r)) 



mez a?^a° 



= E E ll- a) (->--fr)|| LOO{/m) TV; e/m (^)(r,.-^r)) 

A?^A° m6J{j>,e) 

+ E E f T)ll^(7 ro )TK eJro (aW(r,.- Mr)) 

A j VAOmG(J(j,r,£)) c 

<2 ^ ||^')(r,-)||L<-( ( _ 00l _i ltfi)U (j M - i ,oo))2V(ff (fc) (r,-)) 

+ 2 V ||^)(r,-)|| £ =c3V A o_ A o A o„ A o (a«(r,-)) 

A o^o (-^-2+^^.-^+2+^1^) 



< 2nTV>«(0, •)) ei + 2 ^ ||<7 W) ||l<» e x . 

j 

Thus, we obtain that, for any r G (0,To), 

JToE E \\^\r r -^T)\\ L ^ iIm) TV xeIm (a^(T,x-^r))=0. 

Meanwhile, let 

G(r) := C(n) max ||a^(r, maxTV(^')(r, •)) < C(n)||C/ 1 (-)||Lo Ty(C/ 1 (-)), 

J 3 
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where the constant C(n) depends only on the dimension n. Then 

£ £ \\a^\r,.-^T)\\ LX{Im) TV xeIm (a( k HT,x-^r))<G(T). 

Notice that 

/•To 

/ G(r)dr < C^ToWU^-^TV^i-)) < oc. 
J o 

Then, by the dominant convergence theorem, we have 

I r E E - ^r)|| i0 o (/m) TK. e/m (a( fc )(r,x - ^r))dr -4 0. 

U /n c c 

Therefore, we conclude 



lim 

e- 



sup \\U £ (t, ■) -U^(t,-)\\ L i =o(e) whene^O. 

0<t<T /e 



□ 
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